Discussing direct search of dark matter particles in the Minimal Supersymmetric 
extension of the Standard Model with light neutralinosj 
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We examine the status of light neutralinos in an effective Minimal Supersymmetric extension of 
the Standard Model (MSSM) at the electroweak scale which was considered in the past and discussed 
in terms of the available data of direct searches for dark matter (DM) particles. Our reanalysis is 
prompted by new measurements at the Tevatron and B-factories which might potentially provide 
significant constraints on the MSSM model. Here we examine in detail all these new data and show 
that the present published results from the Tevatron and B-factories have only a mild effect on 
the original light neutralino population. This population, which fits quite well the DAMA/LIBRA 
annual modulation data, would also agree with the preliminary results of CDMS, CoGeNT and 
CRESST, should these data, which are at present only hints or excesses of events over the expected 
backgrounds, be interpreted as authentic signals of DM. For the neutralino mass we find a lower 
bound of 7-8 GeV. Our results differ from some recent conclusions by other authors because of a 
few crucial points which we try to single out and elucidate. 



I. INTRODUCTION 

Much interest has recently been raised by some new 
hints of possible signals of dark matter (DM) particles in 
experiments of direct detection (CDMS [l|, CoGeNT 0, 
CRESST which previously reported upper bounds 
only. These hints are in fact only constituted by excesses 
of events over what would be expected from backgrounds. 
What is intriguing is that these events, if actually due 
to DM particles with a coherent interaction with the 
atomic nuclei of the detector material, would be concen- 
trated in a physical region which, expressed in terms of 
the WIMP mass and of the WIMP-nucleon elastic cross- 
section, agrees with the physical region established with a 
high statistical significance by the DAMA Collaboration 
from a measurement of annual modulation over 13 yearly 
cycles with the DAMA/Nal and the DAMA/LIBRA ex- 
periments 3. 

These results have prompted a large number of phe- 
nomenological papers focussed on WIMPs with a light 
mass (around 10 GeV) and a WIMP-nucleon elastic 



cross-section of order (10 



-40 



10 



-41, 



cm , whereas pre- 



vious theoretical and experimental considerations were 
prevalently directed toward physical regions with much 
higher masses and lower cross-sections. Turning to a spe- 
cific candidate, it has now become common to consider 
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neutralinos of light mass (~ 10 GeV). 

Actually, already long ago in Ref. [H] it was stressed 
that, in case of R-parity conservation, a light neutralino 
(i.e. a neutralino with m x < 50 GeV), when it happens 
to be the Lightest Supersymmetric Particle (LSP), con- 
stitutes an extremely interesting candidate for the dark 
matter in the Universe, with direct detection rates acces- 
sible to experiments of present generation. In Ref. [H] it 
was also derived a lower bound of m x ~ 7 GeV from 
the cosmological upper limit on the cold dark matter 
density. The theoretical framework, considered in Ref. 

which allows neutralinos with a mass in the range 
7 GeV < m x < 50 GeV is an effective Minimal Super- 
symmetric extension of the Standard Model (MSSM) at 
the electroweak (EW) scale, where the usual hypothesis 
of gaugino-mass unification at the scale of Grand Uni- 
fication (GUT) of the SUGRA models, is removed; this 
effective MSSM is very manageable, since expressible in 
terms of a limited number of independent parameters. 
This model is the theoretical basis we also adopt for the 
phenomcnological investigations presented in this work 
(for simplicity, the model which entails low-mass (m x < 
50 GeV) neutralino configurations within this effective 
MSSM will hereby be dubbed Light Neutralino Model 
(LNM)); its main features are briefly summarized in Sect. 

eh 

When the DAMA Collaboration published their ex- 
perimental results collected with a Nal detector of 100 
kg over 7 annual cycles in Ref. Q it was proved 
that indeed the population of light neutralinos 0] fitted 
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well these data. The possible interpretation of the an- 
nual modulation results in terms of light neutralinos was 
further confirmed in Refs. d, when the channeling 
effect was taken into account in the experimental analy- 



lOj and the first DAM A/LIBRA combined data were 
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presented 

We recall that the present collection of data by the 
same Collaboration [1] amounts to an exposure of 1.17 
ton x year, with an evidence for an annual modulation 
effect at 8.9 a C.L.. This extended collection of data en- 
tails that the DAM A/LIBRA annual modulation regions 
in the plane WIMP mass - WIMP-nucleon scattering 
cross-section, reported in the figures of Ref. 0] with a 
comparison with the theoretical predictions within LNM, 
essentially maintain their shapes, but with a statistical 
significance increased from 6.5 a to 7.5 a (see Ref. Q for 
a detailed definition of these regions). 

As mentioned above, other experimental collaborations 
have recently reported some events that might be in ex- 
cess of the expected backgrounds. Two candidate events 
for dark matter, which would survive after application of 
various discrimination and subtraction procedures, were 
reported by the CDMS Collaboration fl. In Ref. [l2| 
it is shown that, should these events be actually due to 
dark matter, they would be compatible with light neu- 
tralinos and the DAMA results. Likewise, compatible 
with the LNM and the previously quoted experimental 
data, would be the excess of bulk-like events reported 
by the CoGeNT Collaboration again in case these 
might actually be significant of a DM effect. Most re- 
cently, also the CRESST Collaboration has presented re- 
sults which, if interpreted as due to DM particles, would 
point to WIMPs with a mass < 15 GeV and a scalar 
WIMP-nucleon cross section in the ballpark of the pre- 
vious experiments 0]. 

As previously mentioned, the new experimental results 
of Refs. [H-Q, combined with the DAMA data, have 
recently triggered a large number of phcnomenological 
considerations on coherently-interacting WIMPs of light 
mass (see for instance Refs. 13H24J), with emphasis for 



the mass range 7-10 GeV. We wish to stress that, though 
this is a very interesting mass interval, already pointed 
out in Refs. [1,0^1, the physical region compatible with 
the DAMA results alone and accessible, in particular, to 
relic neutralinos is much wider; its specific size depends 
on a number of astrophysical features (e.g. WIMP distri- 
bution function in the galactic halo) and on the detector 
response (role of the channeling effect in detectors which 
use crystals 25|). For instance, it was shown in Ref. Q 
that for the case of a WIMP halo distribution given by 
a cored-isothermal sphere the extended mass range is 7 
GeV < m x < 60 GeV. In view of the high significance 
of the annual modulation results versus the still prelimi- 



nary character of the hints of experiments in Refs. Q-{^> 
we will consider in the present paper the status of relic 
neutralinos in the whole WIMP light-mass range 7 GeV 
< m x < 50 GeV, as we did in our previous publications 
on light neutralinos. 

Our present investigation is mainly focussed on the 
role that recent measurements at the Tevatron and at 
the B-factories BaBar and Belle can have in provid- 
ing constraints on supersymmetric models more stringent 
than the ones previously considered in Ref. 0). Specifi- 
cally, we consider the new data concerning the following 



processes: a) the decay B s 



/!+, the top-to- 



bottom quark decay with emission of a charged Higgs 
(t — > b + H + ), and the searches for neutral Higgs bosons 
into a tau-lepton pair, at the Tevatron; b) the rare decays 
B ->• t + v T and B ->• D + r + v T (and B -> D + I + v u 
where I = e,/i), at the B-factories BaBar and Belle. 

These measurements have potentially a significant role 
in constraining the supersymmetric parameter space in 
the region of high tan/3 and of light Higgs masses and 
consequently in determining the allowed ranges for a 
number of crucial quantities: neutralino relic abun- 
dance, lower bound on the neutralino mass and elas- 
tic neutralino-nucleon cross-section. Here the nature of 
these constraints is critically analyzed, taking also into 
account the fact that many of them are still affected by 
sizable uncertainties. Their impact in constraining the 
parameter space of our supersymmetric model is investi- 
gated in a twofold way: by analytic investigations and by 
detailed numerical analyses. We believe that the analytic 
derivations are necessary: i) to clarify how the relevant 
physical quantities mentioned above depend on the pa- 
rameters of the LNM, ii) to establish the impact that 
each specific constraint has on these model parameters 
and then consequently on the physical quantities, iii) to 
direct the numerical analyses to the regions of the pa- 
rameter space which are of most relevance, iv) to inter- 
pret the outcomes of the numerical evaluations correctly 
and, finally, v) to serve as a guide for an educated guess 
about how the present situation could evolve, as new ex- 
perimental limits on supersymmetric parameters from ac- 
celerator experiments and other precision measurements 
might become available in the future. 

In the course of our discussion we will also comment 
on how our results compare with some of the outcomes of 
other recent papers where numerical analyses of the pre- 
viously mentioned constraints have been discussed (see, 
for instance, Refs. flj Ej. I2II-I2I W\ ) . 

The scheme of our paper is the following. In Sect. |TT] 
the main features of our effective supersymmetric model 
are presented. In Sect. IIIII we derive from the analytic 
expression of the neutralino relic abundance the lower 
bound for the neutralino mass in a form which displays 
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FIG. 1: Scatter plot of the light neutralino population for the 
LNM-.4 scan, shown in the plane m x — tan /3. The (red) solid 
line represents the analytic bound from the neutralino relic 
abundance given in Eq. ((13} . 
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FIG. 2: Scatter plot for T(Z -> xx) versus m x for the LNM- 
A scan. The (red) solid horizontal line denotes the present 
experimental upper bound to the invisible width of the Z- 
decay into non Standard Model particles. 



its dependence on the main model parameters. Then 
in Sect. [V] we derive an approximate expression for 
the neutralino-nucleon cross-section, which provides an 
easy estimate for this quantity. In Sect. IVII we give an 
overview of the most relevant constraints which can have 
an impact on the ranges of our model parameters. Our 
evaluations are compared to the current results from ex- 
periments of direct DM particles in Sect. I VIII We draw 
our conclusions in ScctllXI 



II. 



LIGHT NEUTRALINOS IN AN EFFECTIVE 
MSSM (LNM) 



The supcrsymmctric scheme we employ in the present 
paper is the one described in Ref. jjj: an ef- 
fective MSSM scheme (effMSSM) at the electroweak 
scale, with the following independent parameters: 
Mi, M 2 , M 3 , //, tan/3, m^, m^, mr and A. We stress that 
the parameters are defined at the EW scale. Notations 
are as follows: Mi, M 2 and M 3 are the U(l), SU(2) and 
SU(3) gaugino masses (these parameters are taken here 
to be positive), [i is the Higgs mixing mass parameter, 
tan/? the ratio of the two Higgs v.e.v.'s, rriA the mass 
of the CP-odd neutral Higgs boson, rrtg is a squark soft- 
mass common to all squarks, mr is a slcpton soft -mass 
common to all sleptons, and A is a common dimcnsionlcss 
trilinear parameter for the third family, A~ h = A^ = Arriq 
and Af = Amj (the trilinear parameters for the other 
families being set equal to zero). In our model, no gaug- 



ino mass unification at a Grand Unified scale is assumed. 

The following experimental constraints arc imposed: 
accelerators data on supcrsymmctric and Higgs boson 
searches at the CERN e+e~ collider LEP2 the up- 
per bound on the invisible width for the decay of the Z— 
boson into non Standard Model particles: T(Z — > xx) < 
3 MeV 29, 30| (the role of this bound will be discussed 
in Sect. IIII Al) : measurements of the b — s- s + 7 de- 
cay process 2.89 < BR(b -> s 7 ) • 10 4 < 4.21 is 
employed here (this interval is larger by 25% with re- 
spect to the experimental determination 31| in order to 
take into account theoretical uncertainties in the super- 
symmetric (SUSY) contributions 32j to the branching 
ratio of the process (for the Standard Model calcula- 
tion, we employ the NNLO results from Ref. [33j)); the 
measurements of the muon anomalous magnetic moment 



2)/2: for the deviation, Aa p 



cxp 



-.the 



fl M — ISM ~ ")!"• u " c <-10via.u11_.il, -auj, ^ 

of the experimental world average from the theoretical 
evaluation within the Standard Model we use here the 

derived from the 
data (the 



(2 a) range 31 < Ao^ • 10 11 < 479 



latest experimental [34| and theoretical 



supcrsymmctric contributions to the muon anomalous 
magnetic moment within the MSSM are evaluated here 
by using the formulae in Ref. [36|; the constraints on 



the SUSY parameters obtained from the searches for the 
neutral Higgs-boson at the Tevatron 37-3^]; the upper 
bound (at 95% C.L.) on the branching ratio for the decay 
B s -> /!+ + n~: BR{B S -> < 5.8 x 10" 

the constraints related to AM b s = Mb — Mi 
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and 
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FIG. 3: Scatter plot of the absolute value of the Wilson co- 
efficients \C S \ for each SUSY contribution to the BR(B S — > 
fi + /j,~) (colored points), compared to the full calculation of 
the dominant term (black points) (for expressions of these 
quantities see Refs. [56H59I ]). for the LNM-„4 scan. The sign 
of each term is indicated in parenthesis as "[+]" and "[—]". 
The values of \C§\ are plotted as functions of the neutralino 
mass m v . 



The role of these two last categories of bounds in con- 
straining the LNM is elucidated in Sect. IVI A[ while the 
constraints from the searches for neutral Higgs-bosons 
at the Tevatron are discusses in Sect. IVI B[ The cosmo- 
logical upper bound on Cold Dark Matter (CDM), which 
is also implemented in our calculations, is discussed in 
Sect. IIIII Other possible constraints from the Tevatron 
and B-factories are analyzed in detail in Sect. IVII 

The linear superposition of bino _B, wino and of 

the two Higgsino states H^, which defines the neu- 
tralino state of lowest mass m x is written here as: 



X = axB + a 2 W {3) + a 3 H° + a 4 H°. 



(1) 



The properties of these states have been investigated in 
detail, analytically and numerically in Ref. [431 ] for the 



case when the smallest mass eigenstate Xi ( or X i n short) 
is light, i.e. m x = m Xl < 50 GeV. Of that analysis we 
report here only the main points which are relevant for 
the present paper. 

We first notice that the lowest value for m x occurs 
when: 



Mr < \n\ , M 2 



(2) 



since the LEP lower limit on the chargino mass (jn x ± > 
100 GeV) sets a lower bound on both and M^: 



FIG. 4: The same as in Fig. [3] except that the Wilson coef- 
ficients \Cg\ are plotted as functions of tan ft. 
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FIG. 5: Scatter plot of the ratio between the full numerical 
calculation of BR(B S — > fi + /j,~) and the approximate expres- 
sion of the dominant term BR^{B a —¥ ^i + /i~) given in Eq. 
(|34[) . for the LNM-.4 scan. The scatter plot is shown as a 
function of the neutralino mass m v . 



|/u|,M 2 > 100 GeV, whereas M x is unbound. Thus, \ 
is mainly a Bino; its mixings with the other interaction 
eigenstates are given by: 
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sin Ow sin j3 



(3) 
(4) 
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FIG. 6: The same as in Fig. [5] except that the scatter plot is 
shown as a function of tan /3. 
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/Ksin/3 



Ml sin /3 + \i cos / 



(5) 



where £1 = m x — Mi and 0^ is the Weinberg angle. 
These expressions readily follow from the general analyt- 
ical formulae given in Ref. 43[ by taking tan/3 > 10, as 
consistent with the scenarios discussed below. 

From the above expressions the following relevant 
property holds: x is mainly a Bino whose mixing with 
H° is non-negligible at small \i. In fact, for the ratio 
|«3 |/|eii I one has: 

~ sin(V sin/3 ^ < 0.43 sin/3, (6) 

kl M 

where in the last step we have taken into account the 
experimental lower bound > 100 GeV. 

It is also useful to explicit the connection between the 
neutralino mass m x and the parameter Mi at small m x . 
From the diagonalization of the neutralino mass matrix 
one finds: 



= Mi — sin 9w Mz ( — - cos j3 sin j3 

ai 01 



Employing Eqs. (jSJ - © wc obtain: 



i Mi - sin^ 6 W ^ ( — + Mi 
\i z \ tan p 



2 fj, 



(7) 



(8) 



Mi (1 - 0.16 Mioo) - i- 1 A* 



inn 



35 s 
tan/3, 



GeV, 



Mi 



(1-0.16 Ml0 o) 



,(9) 



where /iioo denotes /U in units of 100 GeV. 
III. COSMOLOGICAL BOUND AND LOWER 
LIMIT TO THE NEUTRALINO MASS 

The neutralino relic abundance is given by: 

9.9 • 10" 28 cm 3 s- 1 



x f 



9*(x f ) 



1/2 



(o-ann") 



(10) 



where (a ann v) = Xf(a ann v) in t, (er ann v) int being the in- 
tegral from the present temperature up to the freeze- 
out temperature T/ of the thermally averaged product of 
the annihilation cross-section times the relative velocity 
of a pair of neutralinos, Xf is defined as 1/ e m x /Tf 
and g*(xf) denotes the relativistic degrees of freedom of 
the thermodynamic bath at Xf. For (a ann ?;) we will use 
the standard expansion in S and P waves: (<7 an nv) ~ 
a + b/(2xf). Notice that in the LNM no coannihilation 
effects are present in the calculation of the relic abun- 
dance, due to the large mass splitting between the mass 
of the neutralino (m x < 50 GeV) and those of sfermions 
and charginos. 

In our numerical evaluations all relevant contributions 
to the pair annihilation cross-section in the denomina- 
tor of Eq. (|10p are included. However, approximate ex- 
pressions for £l x h 2 can be derived analytically; these will 
prove to be very useful to obtain analytic formulae for 
the lower bound for the neutralino mass. 



A. Scenario .4 

We first analyze the case of small values for the mass of 
the CP-odd neutral Higgs boson m^: 90 GeV < wla < 
150 GeV. In this case the main contribution to (cr ann t;)i n t 
is provided by the A-exchange in the s channel of the 
annihilation process x + X — ^ i> + b (one easily verifies 
that when m x < mj,, (er ann v)i nt entails a relic abundance 
exceeding the cosmological bound) . Thus one obtains [H] : 



n x h z 



4.8 • 10" 



GeV 2 ^( X/ ) 1 / 2 a 2 a|tan 2 /3 



4 [1 



(H) 2 K1 2 



[l-m 2 /m 2 ]i/2 (l + e b Y 



(11) 
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where eh is a quantity which enters in the relationship 
between the b-quark running mass and the correspond- 
ing Yukawa coupling (see Ref. 44j and references quoted 
therein). In deriving this expression, one has taken into 
account that here the following hierarchy holds for the 
coefficients a, of x : 

| ai| > |a 3 | > |a 2 |, |a 4 |, (12) 

as easily derivable from Eqs. 

1/2 

As far as the value of g+(xf) is concerned, we notice 
that for light neutralinos Xf ~ 21 — 22, so that neutrali- 
nos with masses m x ~ 6 - 7 GeV have a freeze-out tem- 
perature Tf ~ Tqcdj where Tqcd is the hadron-quark 
transition temperature of order 100 - 300 McV. For defi- 



niteness, we describe here the hadron-quark transition by 
a step function: if Tqcd is set equal to 300 McV, then 

1 /2 

for m x < 6 GeV one has g+(xf) ' ~ 4, while for heavier 



neutralinos g+(xf) 



1/2 



9. In the approximate ana- 



lytic expressions discussed hereafter, we set Xf /g ir (x-f) 1 ^ 2 
= 21/8 (while in the numerical analysis the actual values 
obtained after solving the Boltzmann equation are used). 

In selecting the physical parameter space for relic neu- 
tralinos, a first fundamental constraint to be applied is 
that the neutralino relic abundance does not exceed the 
observed upper bound for cold dark matter (CDM), i.e. 
Q, x h 2 < (flcDMh 2 )m&x- If we apply this requirement, by 
using Eq. (|11[) we obtain the following lower bound on 
the neutralino mass: 



J 



[1 - ml/ml} 1 / 4 
[1 - (2m x )ym\] 



> 7.4 GeV 



V90 GeV/ 



/ 35 



,90 GeV/ \tanP 
I 



0.12 



tl u, 3 



0.12 



max 



(13) 



Here we have taken as default value for (£lcDMh 2 ) max the 
numerical value which represents the 2er upper bound to 



CDM 



derived from the results of Ref. 



For 



€b we have used a value which is representative of the 
typical range obtained numerically in our model: ej, = 
-0.08. 

Eq. ([T3|) . already derived in Ref. Q, is written here in 
a form which shows more explicitly how the lower limit on 
m x depends on the various model parameters. Notice in 
particular that this lower bound scales (roughly) as m\ 
and (tan/3) -1 . It is obvious that the precise value for 
the lower limit has however to be ascertained by numer- 
ical evaluations which take into account all the intricate 
interferences of the various physical constraints over the 
model parameters. 

The right-hand-side of Eq. (|13p can be expressed com- 
pletely in terms of the independent parameters of our 
SUSY model. In fact, by using Eq. ((4]) at large tan/3 
(sin /3 ~ 1) and taking into account that, because of Eq. 



1 



03, we can rewrite a 2 a 2 as 



1 U 3 



sin 2 9w m 2 z fj, 2 



0.19^ 00 



(M 2 



sin t>w 



m 2 z ) 2 



(M 2 



100 



0.19) 2 



(14) 



From this formula and the LEP lower bound |/^| > 100 
GeV, we obtain (a\a\) max < 0.13. 

An upper limit on a 2 and then on the product a 2 a 2 
is also placed by the upper bound on the width for the 
Z-boson decay into a light neutralino pair. This decay 



idth is given by [46|, |47 1 



T{Z XX) = li^| M|[1 ~ (2-x) 2 /M 2 ] 3 / 2 (a 2 -a 2 ) 2 . 

(15) 

Taking into account that a 2 ^> a 2 : 

r(Z -> XX) = 166 MeV [1 - (2m x ) 2 /A/ 2 ] 3 / 2 of. (16) 

Denoting by T(Z — > xx)ub the upper bound to the invis- 
ible fraction of the Z-decay width, we finally obtain: 



4 < 



154 MeV 



1/2 



(17) 



where we have used for the neutralino mass the value 
Tci-y ~ 10 GeV. 



29 



If we take conservatively T(Z — > xx) < 3 MeV 
30(, from Eq. (fl?)) we find a\a 2 < 0.12, a value which 



is extremely close to the upper bound (afa|) max = 0.13 
derived above from Eq. (|14p and the experimental lower 
limit on |/x|. The value a\a\ = 0.12 is the reference value 
for a 2 a 2 employed in Eq. (fTB"]) . Notice that the upper 
bound on a 2 , placed by the invisible width for Z —> x+X> 
scales with the square root of the upper limit on this 
quantity (see Eq. (jXTJ)) ; furthermore, the lower bound 
on m x scales as (afag) -1 / 2 (see Eq. fTS"]) ). Thus, the 
lower limit on the neutralino mass is only very mildly 
dependent on the actual value of the upper bound on 
T(Z -> xx)- For instance, taking T(Z -> xx) < 2 MeV 
instead of T(Z — > xx) < 3 MeV, would increase the lower 
bound on m x by a mere 10% 
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The properties of the very light neutralinos of cosmo- 
logical interest considered in this Section delineate a spe- 
cific scenario hereby denoted as Scenario A [43|. Its 



main features are strongly determined by the require- 
ment that the neutralino relic abundance satisfies the 
cosmological bound Q x h 2 < (£lcDMh 2 ) max - From the 
approximate formula in Eq. ([lip one finds that: i) rriA 
must be light, 90 GeV < m A < (200 - 300) GeV (90 
GeV being the lower bound from LEP searches) ; ii) tan /3 
has to be large: tan /3 = 20-45, iii) the B — H° mixing 
needs to be sizeable, which in turn implies small values 
of fj,: \n\ ~ (100 - 200) GeV (see Eq. ([13])). As will 
be discussed in Sect. I VI A[ the trilincar coupling is con- 
strained, for neutralinos lighter than 10 GeV, to be in 
the interval \A\ < 0.6, because of the upper bound to 
BR(B S — > /i + + Eq. (fT3")) shows that in Scenario A 
we expect a lower bound on the neutralino mass of the 
order of 7.5 GeV, if the parameter space which defines 
this Scenario is allowed by the bounds on Higgs searches 
and B-physics. We will show in the next Sections that 
this result actually holds. 



B. 



Scenario B 



When niA > (200 — 300) GeV, the cosmological lower 
bound on (cr a nnf)int can be satisfied by a pair annihila- 
tion process which proceeds through an efficient stau- 
cxchangc contribution (in the t, u channels). This re- 
quires that: (i) the stau mass rrif is sufficiently light, 
TOf ~ 90 GeV (notice that the current experimental limit 
is m,f ~ 87 GeV) and (ii) \ is a very pure Bino (i.e. 

I 



(l-a?)~O(10- 2 )). 

The requirement (i) sets a constraint on the quantity 
tan/3, because the experimental lower bounds on the 
sneutrino mass and on the charged slcpton masses of the 
first two families imply a lower bound on the soft slep- 
ton mass: m.- > 115 GeV. Thus, in order to make the 
request rrif ~ 90 GeV compatible with mj > 115 GeV, 
it is necessary that the off-diagonal terms of the slep- 
tonic mass matrix in the eigenstate basis, which are pro- 
portional to [i tan /3, are large. Numerically, one finds 
|u|tan/3 ~ 5000 GeV. On the other side, the condition 
(ii) requires that | £13/01 1 < 10" 1 , i.e., according to Eq. 
0, \a 3 /ai\ ~ sin(Vsin/3(A'/ z /u) < 10 _1 . Combin- 
ing this last expression with the condition |ju|tan/3 ~ 
5000 GeV, one finds that |u| and tan/3 are bounded by: 
> 500 GeV, tan/3 < 10. These bounds are somewhat 
weaker for values of the neutralino mass larger than ~ 
15-18 GeV. 

The previous arguments lead us to introduce Scenario 
B , identified by the following sector of the supersym- 
metric parameter space: Mi ~ 25 GeV, |u| > 500 GeV, 
tan/3 < 10; ?7i,~ > (100-200) GeV, -2.5 < A < +2.5; the 
other supersymmctric parameters are not a priori fixed. 
Within this scenario it follows from Eqs. (|3HS]) that the 
following hierarchy holds for the coefficients af. 



\ai\ > |a 2 |, \a 3 \, |o 4 | 



(18) 



As derived in Ref. [5j the cosmological bound Vt x h 2 < 
(^CDM^Jmai provides the lower bound m x > 22 
GeV, whose scaling law in terms of the stau mass and 
(^cr>M^ 2 )max is approximately given by: 



: [l-ml/m x } 1/4 > 22 GeV ( 



90 GeV 
I 



0.12 



1: 



(19) 



In general, one has conservatively to retain as a lower 
bound to m x the smaller of the two lower limits given 



separately in Eq. (fT3j) and in Eq. (fl9|) . From these equa- 
tions one finds that the lower bound of Eq. (|T3"|) is less 
stringent than the one of Eq. (fT9|) as long as tua < 2mf ■ 
Due to the present experimental bounds on m^tanjS 
and rrif the lower absolute bound is the one derived from 
Eq. (|13[) . We parenthetically note that the lower limits 



m x > (15 - 18) GeV found in Refs. 
the assumption that wia is very large (tua 



are due to 
1 TcV). 



IV. NUMERICAL ANALYSIS OF THE LNM 
PARAMETER SPACE 



In the present paper we are interested in discussing 
neutralinos with very light masses. We will therefore 
concentrate on Scenario A only, and in our numerical 
analyses we perform a scanning of the supersymmctric 
parameter space dedicated to this scenario. We will de- 
note this sector of the LNM parameter space as LNM-A 

The ranges of the MSSM parameters, appropriately 
narrowed in order to explore this scenario are: 10 < 
tan£ < 50, 100 GeV < n < 150 GeV, 5 GeV < M 1 < 
50 GeV, 100 GeV < M 2 < 1000 GeV, 250 GeV < m q < 
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1000 GcV, 100 GeV < mj < 3000 GcV, 90GcV < m A < 
120 GcV, < A< 1. 

In our scenario the trilinear coupling A and the [i pa- 
rameter are both always positive. This is actually due to 
the interplay between the constraints on BR(b — > sj) 
(which requires \iA > 0) and on a M (which requires 
fx > 0). In all the plots shown in the paper, all the ex- 
perimental bounds discussed in Sect. |ll]are applied: i.e. 
invisible Z-width, Higgs searches at LEP and Tevatron, 
BR{b — > S7), muon anomalous magnetic moment a M , 
BR(B S -> /Lt+M"), AM B , S , BR(t bH+), and the up- 
per bound on the cosmological abundance il x h 2 . When- 
ever we refer to "LNM-.A scan", we intend the scan of 
the parameters space defined above, implemented by the 
experimental bounds quoted here. 

As for the scan of the parameter space, we randomly 
sample the above intervals using a logarithmic scale. It 
is worth noticing here that, due to the non-trivial inter- 
play of the different constraints on physical masses and 
couplings, the typical success rate of our sampling for 
obtaining neutralinos with mass less than 10 GeV is of 
order 10 -5 -10 -6 . In particular, in order to populate the 
scatter plots with a sizable number of points and in a uni- 
form way, we have subdivided the above ranges in smaller 
intervals for the neutralino mass m x ~ Mi and run our 
code until a similar number of allowed points were found 
in each sub-range. Less focused scans of the parameter 
space may fail to find allowed configurations, especially 
in the lower range of m x , as seems to be the case with 



some analyses in the literature 21 1. 

All the numerical results shown in the paper refer to 
this special LNM-.A scan. The only exception will be Fig. 
[25l where a more general scan of the effMSSM will be pre- 
sented. In that (unique) case, the parameters will be var- 
ied in the following intervals: 1 < tan j3 < 50, 100 GeV < 
\fi\ < 1000 GcV, 5 GcV < Mi < min(100, 0.5M 2 ) GcV, 
100 GeV < M 2 < 1000 GcV, 80 GcV < m§ < 3000 GcV, 
80 GeV < m t < 3000 GeV, 90 GeV < m A < 1000 GeV, 
-1 < A < 1. The scan of Fig. [25] will therefore in- 
clude both Scenario A and Scenario 23, as well as more 
general scenarios, with heavier neutralinos. Within the 
scan at higher neutralino masses described above we 
have imposed the condition mjvLSP > 1-05 m x (with 
NLSP=sfcrmions,charginos) in order to remove config- 
urations where the relic abundance is determined by 
coannihilations between the neutralino and the Next- 
to-Lightest SUSY Particle (NLSP). Notice that in an 
effective MSSM coannihilations are due to accidental de- 
generacies between uncorrelated parameters. This is at 
variance with the SUGRA scenario, where strong corre- 
lations among the mass of the neutralino and of other 
SUSY particles are expected in particular regions of the 
parameter space. 



Now it is convenient to have a first look at a scatter plot 
for the neutralino population within the LNM-.4. This is 
provided by Fig. [T] where the scatter plot is represented 
in the plane m x — tan /3. In evaluating this scatter plot all 
constraints specifically mentioned in Sect. |TT] (including 
the upper bound on BR(B S — > [j, + fi~) and the cosmo- 
logical upper bound Q x h 2 < (QcDMh 2 ) m ax = 0.12) have 
been applied. Notice that tt x h 2 has been evaluated using 
its full expression, and not simply with its approximate 
version given in Eq. (jllj) . This figure shows how accu- 
rate is the bound given in Eq. (|13[) , which is represented 
by the (red) solid line. From this figure it turns out 
that the lower bound on the neutralino mass is m x > 
7.5 GeV, this value being obtained when tan (3 ~ 40 and 
niA — 90 GeV. We stress that the updated constraint 
on BR(B S — > fi + fi~) induces only a very slight modifi- 
cation in the neutralino mass lower bound of 6-7 GeV 
determined in Ref. j5[- This constraint will be further 
discussed in Sect. IVI Al 

In Fig. [2] we give a scatter plot for T(Z — > \X) ver_ 
sus m x . From this plot one sees that actually the lower 
bound on the neutralino mass changes very little when 
the upper bound on the invisible fraction of the Z-boson 
width, shown by the horizontal solid line, is decreased 
from 3 MeV to about 2 MeV, as previously argued. Some- 
what below this value the impact of the invisible Z-width 
on the lower bound of variation of m x may be substantial. 



V. NEUTRALINO— NUCLEON 
CROSS-SECTION 

In the present paper we analyze the results of present 
experiments searching for direct detection of DM parti- 
cles, under the hypothesis that WIMPs have a dominant 
coherent interaction with the detector nuclei. This is the 
case for neutralinos. Once a specific distribution is as- 
sumed for the WIMPs in the halo, the WIMP-nucleus 
cross-section can be immediately rewritten in terms of 
the WIMP-nucleon cross-section; this is then the cen- 
tral quantity to be analyzed. 

Thus, we turn now to an approximate evaluation of 
the neutralino-nucleon cross-section in the case where 
the interaction process is due to exchange of the lighter 
CP-even neutral Higgs boson h. From the formulae in 
Refs.[5ll. |52| one obtains: 



(nuclcon) 
^scalar 



2^JlM 2 m 2 M 



(20) 



where: 



Fh = (— ai sinOy/ + a,2 cos9w)(0'3 sin a + 04 cos a) 
I h = J2K™ q (N\qq\N). (21) 



9 



The matrix elements (N\qq\N) are meant over the nucle- 
onic state, the angle a rotates H[ ^ and into h and 
H, and the coefficients are given by: 



u — type 
k h 

^d-type 



cos a/ sin/3 , 

— sin a I cos /3 — e d cos(a — (3) tan /? , (22) 



for the up-type and down-type quarks, respectively; e d 
has already been introduced in Eq. f) 1 1 f) for the case of 
the b quark. 

Keeping the dominant terms (couplings of the Higgs 
boson h with the d-type quarks, a ~ tt/2), one has: 



tan /? g d . 



(nuclcon) 
scalar 



6.8 x 10 _7 a 2 a 2 t&n z (3 



A 

4 



(23) 



where: 

g d = [m d (N\dd\N) + m s {N\ss\N) + m b (N\bb\N)}. (24) 
Thus: 



(25) 



I h ~ k'^_ typc [m d (N\dd\N) + m s (N\ss\N) + m b (N\bb\N)} or: 
I 



(nuclcon ) 



scalar 



5.3 x 10 



-41 2 

cm 



"1"3 

0.13 



tan/3V /90 GcV 
35 / V m h 



9d 



290 McV 



(26) 



r 



In this expression we have used as reference value for 
g d the value g d ,vei = 290 MeV employed in our previous 
papers 0, . We recall that this quantity is affected by 
large uncertainties 

(ffd.min/ ' 9d,TefY 



51 



_ with (gd,max/ 9d,xet) = 3.0 and 
0.12 0. Our reference value <?d,rct = 



290 MeV is larger by a factor 1.5 than the central value 
of Ref . |53j , frequently used in literature (see for instance 



Ref. [16|). 



By employing Eq. (|TT]) and Eq. ([25|) we obtain: 



(nucleoli) 
/ ^scalar 



3.3 x 10 



- 39 cm 2 g 2 



2 12 



[1 - (2m x ) 2 /m 2 A 
ml [1 -m 2 /m2]i/2 (1 + e b f 



(27) 



From this expression we find that any neutralino 
configuration, whose relic abundance stays in the cos- 
mological range for CDM ( i.e. (£lcDMh 2 ) m i n < 

n x h 2 < 



1 



with (flcDMh 2 ) 

i: 



0.098 



and ((]coAf/i 2 )max = 0.12) and passes all particle- 
physics constraints, has an elastic neutralino-nucleon 
cross-section of order: 



^ar n) -(2.7-3.4)xl0- 



41 2 

cm 



( 9d V 

V290 McV/ 



[1 - (2m x ) 2 /m 2 



(m x /(10 GeV) 2 [1 - to 2 /™ 2 ] 1 / 2 ■ 



(28) 



A few comments are in order here: 

i) The elastic cross-section cr^""^ 011 ^ is affected by 
large uncertainties because of the uncertainties inherent 
in the effective Higgs-quark coupling constant g d 51 



Actually, ^gcaiar° is subject to an increase by a factor 



of 3.0 or to a decrease by a factor of 8.6 
above; 

(nucleoli) 



as commented 



n J Eq. (|28|) shows that (T^^" U) scales roughly as 
(m x )~ 2 for the range of neutralino masses considered 
here: 
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iii) To establish the range of m x to which Eq. ([28)1 
applies, one simply has to evaluate the lower bound m x 
by using Eq. ([T5J> ; 

iv) The bounds, set by particle-physics measurements, 
on the two parameters tan /? and tua have a strong im- 
pact on the lower bound of m x (see Eq. (flU)) ). but are ei- 
ther uninfluential (in the case of tan /?) or only marginally 
influent (in the case of tua) in the estimate of e^^iar on ^ . 

Furthermore, we wish to notice that also the situation 
when relic neutralinos only provide a fraction of the CDM 
abundance is of great importance. Indeed, as shown in 
Ref. 54j , in a direct detection experiment, relic neutrali- 



nos, whose relic abundance does not saturate the CDM 
abundance (that is, with fl x h 2 < (£lcDMh 2 )min), have 
a response larger than neutralinos of higher relic abun- 
dance. This property is due to the fact that, for sub- 
dominant neutralinos, the direct detection rate has to 
include a factor which appropriately depletes the value 
of the local DM density po when il x h 2 < (QcDMh 2 ) m i n - 
This reseating factor £ = p x /po is conveniently taken 
as £ — min{l 7 il x h 2 /(il C DAih 2 ) m in} [55j ]. Thus, effec- 
tively the relevant quantity to be inserted in the detec- 



tion rate is not simply a, 



nucleoli) 
scalar 



but rather £cr. 



(nucleoli) 
scalar 



. If 



one performs a scanning of the supersymmetric param- 
eter space, it turns out that, at fixed m x , the quantity 

^scaia'i 00 "^' wnen plotted including configurations with a 
subdominant relic density, can provide larger values than 



(nucleoli) 
^scalar 



when the latter is plotted in the case when only 
neutralinos providing the observed DM density are in- 
cluded. This feature is manifest in the numerical results 
presented hereafter. These important rcscaling proper- 
ties are often overlooked in current phenomenological 
analyses of experimental data. 



VI. CONSTRAINTS ON SUPERSYMMETRIC 
PARAMETERS FROM THE FERMILAB 
TEVATRON COLLIDER AND THE 
B-FACTORIES 

Now we discuss some relevant particle-physics mea- 
surements for which there have been sizable improve- 



ments recently, or which might become important in the 
near future as new data become available. We consider 
how each of the present experimental results on searches 
for new physics at the Tcvatron and at the B-factories 
can impact on the LNM-^4 by putting constraints mainly 
on the two crucial parameters tua and tan 0. Once these 
bounds are established, we determine how these limits re- 
flect on the lower bound for the neutralino mass and con- 
sequently on the neutralino-nucleon cross-section. Our 
analysis is performed analytically and numerically indi- 
vidually for each measurement, since the results of the 
various particle-physics experiments do not share the 
same level of reliability; actually, some of them are still 
presented by the experimental collaborations under the 
form of preliminary reports. Thus, in some case it is 
still premature to enforce the corresponding constraints 
at the present stage, though these might possibly become 
relevant in the future. 



A. Search for the rare decay B s 
Tevatron 



at the 



The SUSY contributions to the branching ratio for the 
decay B s — > p + + pT are very sensitive to tan/3, since for 
hig h values of this parameter they behave as tan 6 56- 
59|. Thus, the experimental upper bound on the branch- 
ing ratio for B s — > p + + p~ can potentially put strict con- 
straints on the elastic neutralino-nucleon cross-section 
when this proceeds through a Higgs-exchange, as is the 
case for the cross-section in Eq. (|26p. However, we wish 
to stress here that the actual impact of these constraints 
depends dramatically on the specific SUSY model. 

To clarify this point, we start discussing the features 
of the supersymmetric contributions to the branching ra- 
tio for B s — > p + + p~ which actually go like tan 6 0, 



as derivable for instance from Eqs.(l)-(2) of Ref. [57 1 
(notice however that in the numerical evaluations of 
BR(B S — > p + p~) reported later on, all supersymmet- 
ric contributions are included as given in Ref. [56|). The 
dominant contribution which behaves as tan 6 writes: 



BR<M(B.->n + ir) 



2 12 ^ 3 sin 4 ^ 



C 2 n, 2 

T BM B j B 



m^m t m Xi 
■m 2 w m 2 A 

,2>i 



1 



x\V tb \ 2 \V ts \ 2 tan 6 /3 |Z>(m£//0 - D(m\J p 2 ) 

I 



2 



sin 2 (26»t) x 



(29) 



where tb is the B meson mean life, Mb is its mass, fs 3 is the B s decay constant, m t , mj, m s are the masses of the 
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top, bottom and strange quarks, respectively; Vtb and V ts 
are elements of the Cabibbo-Kobayashi-Maskawa matrix 
and the function D{x) is denned as D(x) = x log(x)/(l — 
x). The structure of BR^(B S -> in Eq. flU} 

is due to the fact that the relevant amplitudes contain a 
one-loop insertion on a quark line, the loop being formed 
by a chargino (of mass m x ± ~ in our models) and a 



stop whose mass eigenvalues are denoted as m 2 (i = 1, 
2); 6i is the rotation angle which comes out when the 
stop squared-mass matrix is diagonalized. If in Eq. ([2"9"]) 
we insert the values: t b = (1.47 ±0.02) x 10~ 12 sec, M B 
= 5.37 GeV, f B , = (210 ± 30) MeV, \V tb \ = 0.88 ± 0.07 
and \V ts \ = (38.7 ± 2.1) x 10~ 3 , we obtain (by using the 
central values for the various quantities): 



BR (6) (B S -> /i+ + ~ 2.70 x 10~ 6 sin 2 (26» t ~) 



V110 GeV 

D{miJ^)-D{mfJ^) 2 



2 /90 GeV\ 4 /tan/? 



35 



(30) 



The uncertainty in the numerical factor in front of the 
RHS of this equation is of about 40%. 

Since the purpose of the present discussion is essen- 
tially illustrative to show in which features the size of 
BR^ in our LNM-„4 may differ from its size in SUGRA 
models, we proceed to some approximations. In the sce- 
nario A of our model (see Sect. [TTJ) is typically small 
(close to the current LEP lower bound \fj,\ > 100 GeV) 
at variance with what occurs in SUGRA models, where 
the SUSY-breaking implies larger values of |/z|; then here 
the ratios m 2 / u 2 are large: m? In 2 3> 1 . 

tl,2 ' tl,2 

Since D(x) ~ — log(a;) when i» 1, one obtains: 



D(ml/n 2 )-D(ml/v 2 ) * log^/mfj .(31) 



Neglecting the contributions of the D-terms, the stop- 
mass eigenvalues are approximately given by: 

m? 2 ^ = m~ + m 2 ± m t (A rriq + ///tan/3). (32) 

In the scenario A, where tan/3 is large and \fi\ is small, 
unless the trilinear coupling is practically null, we have 
|^4| 3> \jj\/(mq tan/3) and consequently: 



2 2 1^1 mt m~ q 



m? + to 2 



(33) 



By inserting Eq. and Eq. ([33]) into Eq. ([3D]) and 

taking into account that in our model sin 2 (20f) ~ 1, we 
obtain: 



BW.^ c U »1H (i^) 2 (l ^) 2 (^) 4 (^) 6 , (34) 



where the numerical coefficient in front of the right-hand- 
side is normalized to the experimental upper bound at 
95% C.L., BR{B S -> /U+/X") < 5.8 x 10" 8 (this is 
the latest published value by the CDF Collaboration; a 
somewhat smaller value is reported in an unpublished 
CDF Public Note 9892 0). 

Thus, neutralino configurations with a trilinear cou- 
pling parameter: 



1 mi 



< 0.36 



(35) 



14 m t rriq 

are compatible with the constraint imposed by the up 



per bound on the branching ratio of the B s 



/' 



process. In the last step we have taken rriq ~ 1 TeV. 
We recall that, because of the uncertainties involved in 
the determination of the numerical factor in Eq. (J3UJ), 
the numerical coefficients in Eq. ([35]) are affected by an 
uncertainty of about 20%. In addition, the approxima- 
tion of Eq. ([3T]) is only partly valid, since in our scenario 
the arguments X{ = m 2 / [i 2 of the D{x) function are 
large but not exceedingly large and a further correction 
is at hand: a careful analysis shows that the value in Eq. 
([341 is actually reduced by a factor rj = (0.75 -t- 1), for 
Xi = (10-7- 100). This implies that the range on |^4| given 
in Eq. (|35p can extend to values larger by a factor rf 1 ! 2 , 
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i.e. \A\ < (0.36-7-0.42). 

We will see now that the situation is even more favor- 
able than the one depicted in Eq. once the role 
of the other SUSY contributions concurring to the full 
calculation of BR(B S —¥ is taken into account. 
In Figs. |3] and @] we show the absolute value of the 
Wilson coefficients for each SUSY contribution to the 
BR(B S — > // + yu~) (colored points), compared to the full 
calculation of the dominant term (black points) (for ex- 
pressions of these quantities sec Rcfs. [56l459j ). for our 
full scan in LNM-A The sign of each term is indi- 
cated in parenthesis as "[+]" and "[—]". We see that 
also some other terms (notably, the W-boson, the Higgs 
and the penguin diagrams) can contribute significantly, 
some with opposite signs, to the total branching ratio. 
The result of the full calculation when compared with 
the approximate expression of Eq. (|34|) is shown in Figs. 
[5] and [51 For models with light neutralinos, especially 
in the case of m x < 10 GeV, the full calculation can be 
smaller that the approximate one by up to a factor of 
three. 

Fig. [7] shows the correlation between tan (3 and the 
trilinear parameter A in the LNM-.4 scan. We remind 
that here, as in all figures, the bound BR(B S /i + /j~) < 
5.8 x 10~ 8 is implemented. We see that, because of the 
competition among various contributions, the range of A 
for light neutralinos (red circles) is wider than the one 
derived in Eq. (J35J) and can extend up to about 0.6. 

Taking into account the requirements on \A\ used 
in deriving the previous analytical approximations, i.e. 
\A\^> \fi\/(mq tan/3), and the upper bound of Eq. (j3jj) , 
we find that \A\ has to satisfy the conditions: 



N. Fomengo, S. Scope!, A. Bottino (SOLO J 



Ml 



tan /3 



« \A\ « 



m t 



(36) 



i.e. a hierarchy which is naturally realized (i.e. no fine- 
tuning is involved) in our model where the values of the 
parameters are defined at the EW scale and not induced 
by SUGRA conditions. 

A demonstration of how this hierarchy is actually real- 
ized is provided by the numerical results in Fig. [8j The 
three separate regions correspond to the numerical val- 
ues for the three quantities specified in the picture and 
in Eq. ((561 . The (coloured) circles denote the neutralino 
configurations with m x < 10 GeV. 



As discussed in the literature [4l|, |42j a supersym 



metric contribution leading to an increase of the de- 



cay rate of the process B s — > /i + + 
to a decrease of the difference AMb, 



fi is correlated 



Ah 



M 



compared to the value expected in the Standard Model, 
AMjf A /- I n Fig- IU we show a scatter plot of the ratio 



Ram b ,s 
the LNM-.4 scan 



AM§ U 3 SY /AM 



SM 
B,s 



as a function of 



for 




FIG. 7: Scatter plot which shows the correlation between 
tan /3 and the trilinear parameter A, for the LNM-.4 scan. 
Black points refer to m x > 10 GeV, red circles to light neu- 



tralinos with m x < 10 GeV 



0.80 ± 0.12 [42J, which at 95% C.L. implies an allowed 
range 0.57 < Ram b , s < 1-03, one sees that the quantity 
AMb, s does not imply any additional constraint on the 
LNM parameter space. 

To summarize the results of the previous discussion, 
we can say that in the LNM the upper bound on the 
branching ratio of the B s — s- + process is not sig- 
nificantly constraining, due to the relatively small values 
of the chargino mass and a small splitting between the 
two stop masses. These are situations which are nat- 
urally obtained in our LNM, contrary to the situation 
that occurs in SUGRA-like models, as the one consid- 
ered in Ref. [15 L It is then erroneous to interpret the 
results of Ref. 15| as valid also for a generic effective 
MSSM with light neutralinos, as is sometimes done in 
recent literature (see for instance Refs. 2ll. |23|). 



B. Search for neutral MSSM Higgs bosons at the 
Tevatron 

The Tevatron is expected to have a good sensitivity to 
the search for Higgs neutral bosons in the regime of small 
niA and large tan /3, since in this region of the super- 
symmetric parameters the couplings of the neutral Higgs 
bosons <f> = h, A, H to the down-fermions are enhanced. 
This has prompted searches by the CDF and DO Col- 
laborations for the neutral Higgs bosons which decay as 
(j) — y bb or (J) — )■ tt (for an updated review see Ref. [61(). 



Taking into account that i?A 



We report here the results of these Collaborations in 



13 



N. Fornengo, S. Scopel, A. Bottino (SOLO J 



N. Fornengo, S. Scopel, A Bottino (2010) 



o 

0) 




0.01 



;(a# 6 / / / / ' / 




•/ /7 ///\/ / 


/ - 


A sr\ / / / /v-- 


J? y /- 


-yf 'a- /\ / w / / 


/ / - 


f/i/ p. \( / / 

i/ / A \/ 12 / (b) / 

/ / / Vv^ 4 


jr y 
/ / 


/ / / / ^ / y 
'! / / / y y 


y 

y 















140 

(GeV) 



FIG. 8: Scatter plot which shows the hierarchy of Eq. (|36[) 
in the LNM-.4 scan. The three separate regions correspond, 
from top to bottom, to the numerical values for the three 
quantities rriq/mt, A and fj,/(iriq tan/3). Coloured circles de- 
note the neutralino configurations with m x < 10 GeV. 
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FIG. 9: Scatter plot of the ratio Ram b , s = AM§ U S SY /AMI, 
as a function of m x , in the LNM-.4 scan. 



terms of upper bounds on tan /3 versus rriA- These bounds 
are displayed in Figs. [TUl[TTl[T2]as piccewise linear paths. 

The DO Collaboration has determined upper bounds 
for the production rate of the process pp_ 



(inclusive t + t~ production) in Rcf. |37| and for the 
t + t~ production in association with a b quark in Rcf. 
[3^ ] and then converted these bounds into upper limits 
for the SUSY parameters tan j3 and vtla- These limits are 



FIG. 10: Upper bounds in the uia - tan f3 plane, derived from 
searches of the neutral Higgs boson at the Tevatron: line (a) is 
from Ref. [13, line (b) from Ref. line (c) from Ref. 
The dot-dashed line (d) represents the preliminary bound 
given in Ref. (kJ. The solid bold lines labeled by numbers 
denote the cosmological bound Q x h 2 < (Qcdm fc 2 )mai for a 
neutralino whose mass is given by the corresponding number 
(in units of GeV), as obtained by Eqs. (jlll I13|) with et, = 
—0.08 and (f2cr>iu ft 2 )max = 0.12. For any given neutralino 
mass, the allowed region is above the corresponding line. 



represented in Figs. [TUl [Tl] [12] by the line (a) (from Ref. 
[33]) and the line (c) (from Ref. [H]). 

The CDF Collaboration has reported upper limits for 
the production rate of the inclusive t + t~ production in 
Ref. 39[ with ensuing upper bounds represented by the 
line denoted as (b) in Figs. Hp ] \U \ [Tl These results 
supersede the stricter bounds found by the same Collab- 
oration in a previous analysis where tan /3 < 40 at = 
90 GeV [fi3]. 

The DO and CDF Collaborations have also presented 
a combined analysis of their searches for Higgs into the 
inclusive t + t~ channel which provides upper bounds on 
tan /3 displaying a sharp variation at small values of m.A'- 
tan /3 < 30 - 31 at m A = 90 GeV and tan f3 < 44 - 46 at 
ra A = 100 GeV (see line denoted as (d) in Figs. [TU| 

mis). 

It is worth remarking that the derivation of the bounds 
on the SUSY parameters from the experimental data re- 
quire the use of a specific supersymmetric model. The 



one employed in Refs. 37H39I. |62j, |63j is different from 
the LNM; in the scenario A of our model, because of the 
typical small values of the parameters \x and m^, the 
bounds on tan (3 and wia might be more relaxed (see ar- 
guments in Sect. 3.2.1 of Ref. 64[). Notice also that 
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FIG. 11: Upper bounds in the mA - tan/3 plane, derived 
from searches of the neutral Higgs boson at the Tevatron: 
line (a) is from Ref. [3?]], line (b) from Ref. [3i|, line (c) from 
Ref. [38j]. The dot-dashed line (d) represents the preliminary 
bound given in Ref. [g^]. The scatter plot refers to the light 
neutralino population of the LNM-.4 scan. Black points stand 
for m x > 10 GeV, while the red circles for m x < 10 GeV. 



the results of Ref. 63[ are still presented as an (unpub- 



lished) preliminary report. Thus in our analysis we only 



employ the bounds of Refs. j37H39| , which taken together 



disallow the region depicted in yellow in Figs. [T01 ITTl [T21 
In Fig. [10] we display (as continuous curves in black) 
the lines where £l x h 2 , calculated with Eq. (fTTj) , is equal 
to (f^C-DAf ft- 2 )max = 0.12 at the fixed value of m x indi- 
cated (in units of GeV) along each curve; for the other 
parameters the following values are used: Xf /g ir (xf) 1 ^ 2 
= 2.63, a\a 2 = 0.12, e b = —0.08. For a given value of 
m x (masses from 6 GeV to 20 GeV are considered here) 
the region below the relevant curve is disallowed by the 
cosmological upper bound on QcDAih 2 . By comparing 
the continuous (black) curves with the Tevatron limits, 
one sees what is the impact of these limits over the al- 
lowed range for the neutralino masses. In particular, one 
notices that the (yellow) forbidden region in compatible 
with neutralino masses down to 7 GeV. Should one in- 
clude the upper bounds of Ref. [(33[ , the lower limit on 
m x would be increased only very slightly to the value of 
about 7.5 GeV. In Fig. [TTJwe instead display the features 
of the scatter plot for the light neutralino population of 
LNM-^4, in the plane ttla — tan/3, once the bounds pre- 
viously discussed are applied. The points denoted by 
red circles refer to configurations with neutralino masses 
lighter than 10 GeV. 

To complete the analysis of the previous Section on 



FIG. 12: The same as in Fig. 1101 except that curves labeled 
by numbers denote the upper bounds from the approximate 
analytic expression of BR^ e '(B s — » /x + /i~) in Eq. (|34[) . for 
m x + =110 GeV. The different lines refer to values of the 
stop-masses splitting parameter 8 — Amtrriq/(m^ + rrit) in 
Eq. (|33|) as given by the corresponding number reported close 
to the lines. For any given value of S, the allowed region is 
below the corresponding line. 



the bounds coming from the B s — > /i + + fj," decay, when 
combined with the Higgs searches at the Tevatron, in 
Fig. [T^] we show the upper bounds obtained by using 
the approximate expression of Eq. (|3~4"|) for BR^\B S — > 
(i + fJ,~). The solid lines refer to the bounds obtained for 
fixed values of the stop-masses splitting parameter S = 
\A\m t mql '(m? + m 2 ), from 5 = 0.02 to 6 — 0.8. For any 
value of J, the allowed region is below the corresponding 
curve. We notice that BR^'{B S — > [i + n~) does not set 
significant bounds as long as i5 is sufficiently small, which 
in turn occurs for small values of \A\. In our scan for 
Scenario A, the values of S naturally range between 0.01 
to 0.6, for configurations with light neutralinos. 



C. Search for charged Higgs bosons in top quark 
decay at the Tevatron 

Supersymmctric models which contain light neutrali- 
nos automatically involve also light charged Higgs bosons 
since, at tree-level, the following relation holds 



2 



/"/? 



b + H+ 



A + m W 



This would make the decay t 



651 ] possible in our LNM. 
A search for the decay t — > b + H + , 



conducted at the 
Tevatron, led the CDF Collaboration [66j to establish an 
upper bound on tan f3 which is a monotonically increasing 
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FIG. 13: Scatter plot of the branching ratio BR(t -> bH + ) 
as a function of the charged Higgs mass m H +, in the LNM- 
A scan. The solid line and the yellow band represent the 
experimental upper bound and its quoted uncertainty [r36| . 
Black points stand for m x > 10 GeV, while the red circles for 
m x < 10 GeV. 
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FIG. 14: Scatter plot of the branching ratio BR(t bH + ) as 
a function of tan j3, in the LNM-.4 scan. Black points stand 
for m x > 10 GeV, while the red circles for m x < 10 GeV. 

function of m#±. In particular, at vn H ± = 120 GeV (i.e. 
at rriA — 90 GeV) this constraint corresponds to tan j3 < 
45-50. The bound on the branching ratio BR(t — > bH + ), 
with its quoted uncertainty (yellow band), is shown in 
Fig. Q2] as a solid line, together with the scatter plot of 
configurations of LNM-^4. The yellow band denotes the 
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FIG. 15: Scatter plot of the branching ratio BR(t — > bH + ) 
as a function of m x , in the LNM-.4 scan. 

quoted uncertainty on the bound [63|. We see that the 
current bounds on the decay t — > b + H + do not impose 
additional constraints on LNM-.4. Figs. H4l and [T5l show 
the correlation of BR(t — > bH + ) with tan (3 and with m x , 
respectively. 



D. B —¥ tv and B ^ D + t + v at Belle and BaBar 

The measurements of the B-meson decays B — > t + v 
and B — > D + r + v (and B — > D + / + i/, where I = e,fx) 
are potentially a way to investigate allowed ranges for 
the two parameters tan j3 and to^ ± . 

However, it is to be noted that at present the uncer- 
tainties affecting the theoretical estimates as well as the 
experimental determinations concerning the class of B- 
meson decays mentioned above imposes a very cautious 
attitude in applying tout court the entailing constraints 
on the SUSY parameters. The situation might evolve fa- 
vorably in the future and thus provide either more solid 
constraints or hopefully a substantial indication of new 
physics. Thus, we devote this section to an analysis of 
these processes more in view of possible prospects for the 
future than for an actual implementation at the present 
stage. 

As for the first process, a convenient quantity to be 
studied is the ratio of the total (SM contribution plus 
extra contributions) branching ratio BR to t{B — > tv) to 
the branching ratio due only to SM, BRsm{B — > tv): 
Rbtv = BR tot (B —> tv)/BRsm(B — s- tv). If, as extra 
contributions, only SUSY contributions are taken, one 
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FIG. 16: Scatter plot of the quantity Rbtv, calculated accord- 
ing to Eqs. l[37 ) -l[38 j) . as a function of tan/3 in the LNM-^4 
scan. Black points stand for m x > 10 GeV, while the red cir- 
cles for m x < 10 GeV. The green horizontal band represents 
the range of Eq. (|4Tj) . 



FIG. 18: Scatter plot of the quantity R(D), calculated ac- 
cording to Eq. (9) of Ref. [7(| as a function of tan ft in the 
LNM-.4 scan. Black points stand for m x > 10 GeV, while 
the red circles for m x < 10 GeV. The green horizontal band 
represents the bottom part of the range of Eq. (|42|) . 
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FIG. 17: Scatter plot of the quantity Rbtv, calculated ac- 
cording to Eqs. (|37|l (|38|) . as a function of m x in the LNM-.4 
scan. The green horizontal band represents the range of Eq. 

(ED- 



FIG. 19: Scatter plot of the quantity R(D), calculated accord- 
ing to Eq. (9) of Ref. [7r3 | as a function of m x in the LNM-.4 
scan. The green horizontal band represents the bottom part 
of the range of Eq. (|42jl . 



finds 67 



where: 



Rbtv — 1 



rrib m T 



C 



NP 



(37) 



tan 2 /? 



NP 



in 



H = 



1 — co ta- n P 



(38) 
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FIG. 20: Rbtv and R(D) as functions of C^ P , for configura- 
tions of the LNM-„4 scan. The curve with the parabolic shape 
and the almost straight line represent the values of Rbtv and 
R(D), respectively, for light neutralinos. The red points de- 
note the neutralino configurations with m x < 10 GeV. The 
two ranges along the vertical axes denote the interval of Eq. 
(141 p for Rbtv (in green, on the left) and the interval of Eq. 
(|42|l for R{D) (in blue, on the right). The two vertical bands 
in yellow denote the ranges of C T NP where the experimental 
intervals of Rbtv and R(D) have a common solution. 



The quantity eo is defined by: 
2 a 



eo 



3 7T 



M 3 nC (ml,ml,M*), (39) 



where: 



Co{x,y,z) 



xy \og(x/ y) + yz log(j//z) + xz log(z/ x) 



{x-y){y- z){z-x) 



For Rbtv wc use here the 95% C.L. range: 
0.44 < R Btu < 2.67, 



(40) 



(41) 



based on the experimental world average BR oxp (B + —> 
t+v) = fl.72±g;ti) x 10" 4 deduced in Ref. [g| from the 
Belle [70( and BaBar 71 1 data and the SM evaluation 
BR SM (B + -> t+v) = (1.10 ±0.29) x 10- 4 (z3, Iz3| (this 
determination for BRsm(B+ — > t+v) is taken conserva- 
tively; estimates by other authors 6t| 74 give slightly 
lower values). 

Fig. [TB] and Fig. [T7] display how the band of Eq. 
(|41[) compares with the population of light neutralinos 
in terms of tan /3 and of m x , respectively. These figures 
show that the present range of Rbtv favors large values 
of tan ft and would not have any impact on the lower 
bound of 7.5 GeV previously established. 



FIG. 21: Scatter plot of CJ^ P as a function of the neutralino 
mass m x for the LNM-.4 scan. The two horizontal bands 
in yellow denote the ranges of C T NP where the experimental 
intervals of Rbtv and R(D) have a common solution. Blue 
points stand for cosmologically subdominant neutralinos (i.e. 
n x h 2 < 0.098), while the red crosses refer to dominant con- 
figurations (i.e. 0.098 < n x h 2 < 0.122). 



As for the semileptonic decays B — > D + I + is, it 
is convenient to consider the ratio R(D) = BR(B — > 
Dtv) I BR(B — > Dev), since in this ratio many exper- 
imental systematic uncertainties cancel, either partially 
or completely [75| . Also some theoretical uncertainties 
cancel out in this ratio [7f| , 

An experimental determination for R(D) is provided 
by the BaBar Collaboration: R(D) = (41.6±11.7(stat)± 
5.2(syst)) x 10" 2 [zi|. A determination by the Belle 



77 



Collaboration is given in the report of Ref. 
R(D) = (60 ± 14(stat) ± 8(syst)) x 10" 2 . However, this 
value has been obtained by using the Belle determina- 



-> D°t+v) = (1.51 



^o'39( s t a t) 



+0.24 
19 



(syst) ± 



tion BR(B+ 

0.15) x 10~ 2 , which has recently been superseded by 
the new Belle determination BR(B + — > D°t + v) = 
(0.77±0.22(stat)±0.12(syst)) x 10~ 2 [zi|, smaller than 
the previous one by a factor of two. In view of the lack of 
an updated value for R(D) provided by the Belle Collab- 
oration and of the ambiguities which may rise in treating 
the uncertainties in the various branching ratios concur- 
ring in the determination of R(D), in our discussion we 
only take into account the BaBar value previously men- 
tioned. Notice however that the new Belle value for the 
branching ration of B + — > D°t + v approaches now con- 
siderably the BaBar value for the branching ratio of the 
process B + D°t + i> fr^ . 

From the BaBar determination R(D) = (41.6 ± 
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FIG. 22: Scatter plot of CJ^ P as a function of tan /3 for the 
LNM-„4 scan. The two horizontal bands in yellow denote 
the ranges of CJ^p where the experimental intervals of Rbtv 
and R(D) have a common solution. Black points stand for 
m x > 10 GeV, while the red circles for m x < 10 GeV. 



11.7(stat) ± 5.2(syst)) x 10~ 2 we obtain the 95% C.L. 
range: 



13.5 x 1(T 2 < R(D) < 69.7 x 10" 



(42) 



Fig. [18] and Fig. [19] display the scatter plot for the 
quantity R(D) evaluated by using the Eq. (9) of Rcf. 
|76j . versus tan/3 and m x , respectively. We see that the 
experimental range of R(D) given in Eq. (|42p is compat- 
ible with the light neutralino configurations. 

In Fig. [50] we display Rbtv and R(D) as functions of 
the quantity CJf P defined in Eq. (|38|) . The curve with 
the parabolic shape represents Rbtv as given by Eq. ([57]) 
and the almost straight line gives R(D), calculated with 
Eq. (9) of Ref. [zf|. The parts colored in red pertain to 
neutralino configurations with m x < 10 GeV. 

Notice that, should the indication of Fig. [2U] be 
taken strictly, one would deduce for the quantity Cj^p 
two ranges of compatibility: —0.51 < C T NP < —0.32 
and -0.041 < C T NP < -0.065. From this, by us- 
ing Eq. (|38[) . one finds that the range at large tan/3 
for the quantity tan/3/(m#±/120 GeV) is: 30 < 
tan f3/(m H ± /120 GeV) < 40. 

This in turn would have some impact on the features of 
the neutralino population which substantially contribute 
to the DM abundance (i.e. for which {QcDMh 2 ) m - m < 
Cl x h 2 < (ric_DM^ 2 )max)- This can be appreciated in Figs. 
|2"T1 and [22] where Cn p is plotted against the neutralino 
mass (Fig. [2T) or tan/3 (Fig. [25]). In Fig. [21] the points 



denoted by red crosses refer to cosmologically dominant 
neutralinos, and they are compatible with the preferred 
ranges of Cnp for neutralino masses below 13 GeV (and 
large tan/3) and in the interval 18 — 25 GeV (and low 
tan/3). In Fig. [22] the red circles refer to configurations 
with m x < 10 GeV. Notice that if we take at face value 
the current bounds on the B — > r decays, light neutrali- 
nos would actually be favored. 

Nevertheless, we recall that, for the reasons mentioned 
at the beginning of the present Section, it seems prema- 
ture to enforce these constraints rigidly at the present 
time (see Ref. [69( for similar comments). 

We can conclude the present Section on constraints on 
supersymmctric parameters from the Tevatron collider 
and the B-factories with the following remarks: i) the 
upper bound on the branching ratio for the decay B s — > 
fi + + fx" determined at the Tevatron has a mild effect in 
constraining the LNM, at variance with what occurs in 
SUGRA models, ii) the bounds which are derived from 
the Higgs bosons searches at the Tevatron do not modify 
the previously mentioned lower bound m x > 7.5 GeV, 
iii) the measurements of the rare B-meson decays at B— 
factories have still to be taken with much caution: com- 
bining the present data on B — > r + v and B — > D + t + u 
one derives a range of tan/3/(m#±/120 GeV) which at 
present might only have some effect on the light neu- 
tralino population for m x ~ 15-20 GeV, without modi- 
fying the lower bound on the neutralino mass. 

We finally note that the strict bounds obtained in Ref. 
[lfi | on the WIMP-nucleon scattering cross-section are 
mainly due to: a) the use of constraints on tan /3 de- 
rived from DO-Collaboration [7^] and CDF [(32| data 
which were subsequently superseded (and relaxed) by 



Refs. [37J, [39(, respectively; b) the implementation of 



bounds rcstrictively derived from B-meson decays data 
which suffer from sizable uncertainties and which have 
been in part superseded by the more recent measure- 
ments quoted above. 



VII. CONFRONTING RESULTS FROM 
EXPERIMENTS OF DIRECT DETECTION OF 
DM PARTICLES 

As mentioned in Sect. [V] the neutralino interacts 
mainly by a coherent process with the target nuclei, thus 
the neutralino-nucleus cross-section is conveniently ex- 
pressed in terms of the neutralino-nucleon cross-section 
CT scaia'r and then the relevant quantity to be considered 



is this cross-section multiplied by the rescaling factor £ 
defined in Sect. EJ 

In Fig. [23]wc show the scatter plot representing the su- 
persymmctric configurations of LNM- A and subjected to 
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FIG. 23: c a ( nuclcon ) as a function of the neutralino mass for 

3 SCtLltLr 

the LNM-.4 scan and for gd,icf = 290 MeV (constraint from 
Rbtv not included). The (red) crosses denote configurations 
with a neutralino relic abundance which matches the WMAP 
cold dark matter amount (0.098 < Q x h 2 < 0.122), while the 
(blue) dots refer to configurations where the neutralino is sub- 
dominant (Q, x h 2 < 0.098). The blue-band flag-like region 
denotes the extension of the scatter plot upwards and down- 
wards, when the hadronic uncertainties are included. The 
green shaded regions denote the DAM A/LIBRA annual mod- 
ulation regions [8lJ; the region delimited by the dashed (solid) 
line refers to the case where the channeling effect is (is not) 
included. The two regions are denoted by letters Ai and A2, 
respectively. The yellow hatched regions denoted by letter B 
display the regions (at 68% and 85% C.L.) related to the two 
CDMS candidates [12|. The pink small (horizontally shaded) 
region denoted by letter C refers to the CoGeNT excess of 
events whereas the black straight dot-dashed lines de- 
noted by letter D show schematically a region linked to the 
excess reported by CRESST ||. 



the constraints discussed in Sects. [TTI and [VTl The cross- 
r (nucioon) . g ca i cu i a t e( j w jth its complete expres- 



scction a. 



calar 



sion given in Ref. [52( at a fixed reference set of values 



for the hadronic quantities involved in the neutralino- 
nucleon cross-sections Q (the dominant coupling is 
put at the value ffd.ref = 290 MeV mentioned in Sect. 
|V]) . The (red) crosses denote configurations with a neu- 
tralino relic abundance which matches the WMAP cold 
dark matter amount (0.098 < VL x h 2 < 0.122), while the 
(blue) dots refer to configurations where the neutralino 
is subdominant (tt x h 2 < 0.098). The region covered by a 
(blue) slant hatching denotes the extension of the scatter 
plot upwards and downwards, when the hadronic uncer- 
tainties extensively discussed in Ref. 



51 



are included. 



FIG. 24: The same as in Fig. [23] except that here the con- 
straint from Rbtv is included. 
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FIG. 25: The same as in Fig. [23] except that here the MSSM 
parameter space is scanned beyond the LNM-.4 intervals, as 
specified in Sect. IIV1 in order to include neutralino configu- 
rations of higher mass. 



Notice that the values displayed by the scatter plot 



validate the approximate expression in Eq. (|28[) . They 
differ significantly from the estimates given elsewhere 15, 
16| for the reasons discussed in the previous section. 

To report in Fig. [23] also the results of present ex- 
periments searching for direct detection of DM parti- 
cles, we have to assume a specific model for the distri- 
bution function (DF) of the WIMPs (i.e. the neutrali- 
nos, in our case) in the galactic halo. Among the var- 
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ious possible DFs |80|, we employ here as reference DF 
the one described by the density profile of the cored- 
isothermal sphere (denoted as Evans logarithmic model, 
or Al model, in Ref. (g^) which is given by: 



heavier WIMPs (mass > 20 GeV) was also presented by 
the KIMS Collaboration " 



p(r) 



3R~ 



4ttG (i? 2 + r 2 ) 2 



(43) 



where G is the Newton's constant, u is the local value 
of the rotational velocity and R c is the core radius. For 
R c we use the value R c = 5 kpc. For the parameter vq 
we take the value vo = 220 km sec -1 and for the escape 
velocity the value v esc = 650 km sec -1 . We set p = 0.34 
GeV cm -3 for the total local DM density. 

The green shaded regions denote the DAMA/LIBRA 
annual modulation regions, under the hypothesis that the 
effect is due to a WIMP with a coherent interaction with 
nuclei; the region delimited by the solid line refers to the 
case where the channeling effect is not included, the one 
with a dashed contour to the case where the channel- 



ing effect is included 25 1. These regions represent the 
domains where the likelihood-function values differ more 
than 7.5 a from the null hypothesis (absence of modula- 
tion); they are derived by the DAM A Collaboration 81 1 
from their data referring to an exposure of 1.17 ton x 
year, with an evidence for an annual modulation effect 
at 8.9 a C.L. 0]. 

As mentioned in the Introduction, recently a num- 
ber of experimental Collaborations have reported new 
data consisting of excesses of events (over the expected 
backgrounds) which might represent hints for very light 
DM candidates: CDMS fl, CoGeNT 0, CRESST H. 
Other experimental investigations (the XENON10 [82j 
and XENON100 ^ experiments and the CDMS re- 
analyzes of previously collected data [84|) have led these 



groups to present upper bounds on £,cr^^ on ^ in the same 
range of WIMP masses (around 10 GeV). It is worth not- 
ing that none of these experiments is sensitive to a spe- 
cific signature of DM particles such as the annual modula- 
tion or the directionality. Thus, their detection technique 
and data analysis is based on rather intricate discrimina- 
tion criteria and sizable subtractions, which become more 
and more critical as the analyses of data are extrapolated 
into the range of very low recoil energies. This calls for 
a very cautious attitude both towards a claim of a possi- 
ble signature, in case of excesses of events over expected 
backgrounds, as well as in implementing upper bounds 
rather fragile at the present stage. It is worth stressing 
that a major critical point consists in a reliable determi- 
nation of the actual experimental efficiencies at very low 
recoil en ergies, a problem which is the subject of much de- 
bate 3, 85 - 87 1 . It is beyond the purpose of the present 
paper to enter into these experimental points. We recall 
that an upper bound on ^cr^""^ 001 ^ concerning somewhat 



In Fig. [23] we also illustrate where the regions of the 
possible hints for DM discussed in Refs. [l|-|3| are lo- 
cated. The yellow hatched regions are related to the two 
CDMS candidate events [l|, as derived in [l2| under the 
hypothesis that these events might be due to DM . Also 
displayed in Fig. ESI are the CoGeNT region singled out 
by this Collaboration as due to an excess of bulk-like 
events 0] and a CRESST region which we denote with 
two black dot-dashed straight lines ; this region is meant 
to represent approximately the data reported in Ref. Q 
about 32 signals versus a background estimate of 8.7±1.4, 
compatible with WIMPS with a mass < 15 GeV and a 
WIMP-nucleon cross-section of a few xl0~ 41 cm 2 . 

In Fig. EUwe display the scatter plot of the neutralino 
configurations, when the constraint of Eq. ([4Tj) on Rbtv 
is included. We see that adding this constraint would 
have some impact in depriving the scatter plot of some 
neutralino configurations of significant relic abundance 
in the range 13 GeV < m x < 18 GeV and for m x > 
25 GeV. We remind that in Fig. EH EH only neutralino 
configurations of the special scan LNM-.A are displayed. 

Finally, in Fig. [25]wc show the same scatter plot of Fig. 
1231 where now the MSSM parameter space is scanned 
beyond the LNM-.A intervals, as specified in Sect. IIV[ in 
order to include neutralino configurations of higher mass. 

In conclusion, from the features displayed in Figs. I23[ 
[24l and [25l we derive that: 

i) the light neutralino population agrees with the 
DAMA/LIBRA annual modulation data over a wide 
range of light neutralinos: 7-8 GeV < m x < 50 GeV 

ii) this population is also in agreement with the data 
of CDMS, CoGeNT and CRESST, should these results 
be significant of real DM signals; under these circum- 
stances the range of the neutralino mass would be more 
restricted: 7-8 GeV < m x < (10-15) GeV. 

iii) our results contradict the argument, sometime put 
forward (see for instance Refs. [U, [3]), according to 
which the results of Ref. [3] , which are obtained in the 
context of a SUGRA scenario, imply that the BR(B S — s- 
p + p~) constraint prevents cr^uT^ n ' om reaching the 
size required to interpret current direct-detection data 
in terms of light neutralinos within the MSSM. 



VIII. LINKS TO INDIRECT SEARCHES FOR 
DARK MATTER PARTICLES 

Though indirect searches for relic particles is out of the 
scope of the present paper, some considerations referring 
to indirect effects due to light relic neutralinos are in 
order here. Rather than trying to be complete on these 



21 



topics, we just recall briefly some of the most interesting 
aspects. 

It is known that light relic particles, through their self- 
annihilation processes, can produce a flux of cosmic an- 
tiprotons in excess of the measured one. Indeed, the 
experimental antiproton spectrum is fitted well by the 
secondary component from cosmic rays spallation, calcu- 
lated with the set of the diffusion parameters which is 
derived from the analysis of the boron-to-carbon ratio 
(B/C) component of cosmic rays {8{J, with an estimated 
uncertainty of about 20%. This feature makes the cos- 
mic antiproton flux a potential stringent constraint for 
any exotic astrophysical source of primary antiprotons. 

In Ref.Q the low-mass neutralino populations ex- 
tracted from the DAMA/LIBRA data (depending of the 
size of the channeling effect and on the parameters of 
the halo distribution function) were analyzed in terms of 
the expected effects on the cosmic antiprotons. It was 
concluded that many of these populations arc fully com- 
patible with the current bounds on cosmic antiprotons, 
especially for values of the local dark matter density po 
and local rotational velocity vq in the low side of their 
physical ranges, and for values of the diffusion parame- 
ters of the two-zone propagation model not too close to 
the values of their maximal set [9oj |. A similar analysis 
in the case of the CoGeNT data has been performed in 
Ref.[£H in the phenomenological framework of effective 
DM-quark interactions. 

At variance with cosmic antiprotons, where primary 
and secondary fluxes have very similar behaviours at low 
energies, and can then hardly be disentangled from each 
other, measurements of cosmic antideuterons could pro- 
vide an evidence of light DM particles 92 - 94 1 . In Ref. Q 
it was shown that a sizable number of neutralino con- 
figurations compatible with the annual modulation data 
can generate signals accessible to antideuteron searches 
planned for the next years. 

Also the possibility of investigating light WIMPs at 
neutrino telescop_es has been the subject of specific in- 
vestigations 
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- 99(| . In Rcf.[99j a detailed analysis of 
the neutrino-induced muon signal coming from light- 
neutralino pair-annihilations inside the Sun and the 
Earth is performed and it is shown that, under favorable 
conditions, a combination of the WIMP direct detection 
data and the measurements at neutrino telescopes with 
a low threshold energy could help in pinning down the 
features of the DM particles. 

Other possible signals due to relic light neutralinos are 
mentioned in Ref.[9j. 



IX. CONCLUSIONS 

In this paper we have analyzed the properties of a pop- 
ulation of light neutralinos in an effective MSSM at the 
electroweak scale, already discussed in the past 0, 0], 
in light of new measurements at the Tevatron and B- 
factories which could potentially provide significant con- 
straints in some relevant supcrsymmetric parameters. 

Particular attention is devoted to the branching ratio 
of the process B s — > /i + + /i~ whose experimental upper 
bound entails rather strict constraints on SUGRA mod- 
els. In the present analysis it is shown analitically and 
numerically why this experimental limit has only a mild 
effect on our light neutralino model. 

The light neutralino population, while satisfying the 
cosmological upper bound on cold dark matter, entails 
also a neutralino-nucleon cross-section of the correct size 
to interpret the current experimental results of exper- 
iments for direct detection of dark matter particles in 
terms of MSSM neutralinos. This population, while fit- 
ting quite well the DAMA/LIBRA annual modulation 
data, would also agree with the preliminary results of 
CDMS, CoGeNT and CRESST, should these data, which 
are at present only hints or excesses of events over the ex- 
pected backgrounds, be interpreted as authentic signals 
of DM. For the neutralino mass we find a lower bound of 
7-8 GeV. We have also discussed in detail by how much 
this lower limit would be affected, as more refined and 
solid constraints from the searches on Higgs bosons and 
rare B-decays at the Tevatron and B-factories might be 
derived. It is obvious that great expectations are on the 
outcomes that hopefully will come out from the Large 
Hadron Collider at CERN to bring light to properties re- 
lated to supersymmetry (in Ref. |43[ the perspectives of 
searching for light neutralino of cosmological interest at 
LHC are investigated). 

Our results differ from some recent conclusions by 
other authors; in the course of the presentation of our 
results we have tried to single out and elucidate the main 
points at the origin of these variances. 
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